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ABSTRACT 

Aims. We investigate the morphological, dynamical, and evolutionary properties of the internetwork and network fine structure of the 
quiet sun at disk centre. 

Methods. The analysis is based on a ~6h time sequence of narrow-band filtergrams centred on the inner- wing Can K2v reversal at 
393.3 nm. To examine the temporal evolution of network and internetwork areas separately we employ a double- Gaussian decomposi- 
tion of the mean intensity distribution. An autocorrelation analysis is performed to determine the respective characteristic time scales. 
In order to analyse statistical properties of the fine structure we apply image segmentation techniques. 

Results. The results for the internetwork are related to predictions derived from numerical simulations of the quiet sun. The average 
evolutionary time scale of the internetwork in our observations is 52 sec. Internetwork grains show a tendency to appear on a mesh- 
like pattern with a mean cell size of ~4-5 arcsec. Based on this size and the spatial organisation of the mesh we speculate that this 
pattern is related to the existence of photospheric downdrafts as predicted by convection simulations. The image segmentation shows 
that typical sizes of both network and internetwork grains are in the order of 1.6 arcs. 
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1. Introduction 

The solar chromosphere is the highly structured and very dy- 
namic atmospheric link between the photosphere and the thin 
transition layer to the hot corona. Understanding its physical 
structure is most likely to shed light on outer- atmosphere heating 
mechanisms and the acceleration and composition of the solar 
wind. Well known chromospheric diagnostics are the cores and 
inner wings of strong ab sorption lines lik e Ha a n d Ca ii H & K . 
Sin ce the e arly w ork of I Jensen & Orralll (Il963h . lOrralll (Il966h . 
and ICramI (Il978h decades ago, the resonance lines have been 
extensively employed to explore chromospheric dynamics and 
chrom ospheric structure. For revi ews and refere nces therein we 
refer to lJudge & Pe"te3 (Il998h and lRuttenI (Il999h . 

In low- spatial resolution filtergrams the internetwork re- 
gions appear rather homogeneous and dark, while the highly 
structured and inhomogeneous face of the chromosphere in 
the internetwork is only revealed in high-resolution filtergrams. 
Internetwork areas are intermittently filled with brightenings. 
These internetwork brightenings are most conspicuous in the 
emission feature in the blue wing of the Can K and H lines. 
Thus, they have been named Can K(H)2v grains or simply in- 
ternetwork grains. The investigation of the properties of these 
grains is a diagnostic tool to probe chromospheric conditions in 
internetwork regions. For a detailed description of the K2v grain 
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phenomenon we refer to lRutten & Uitenbroekl (1 19911) and refer- 
ences therein. 

Where does the intermittent chromospheric structure in 
space and time originate from? Observational evidence makes 
clear that the dynamic behaviour and radiative signature of 
the chromospher e is determined by its re sponse to up ward 



travelling waves (IWikst^l et al.L [2000: Judg e et al.L l200lh ex 



cited by the effects of convection (lHoekzema& Ruttenll 19981: 
ISkartlien etaP. l2000l: iHoekzema et al.L l2002h . While the waves 
propagate upwards they steepen into shocks and form high- 
frequency shock sequences pervading the internetwork regions. 
The formation of shocks is strongly suppor te d by t he numeri- 
cal simulations of ICarlsson & SteinI (Il992allbl. Il997h . However, 
because of the one-dimensionality of these simulations, interfer- 
ence with slanted waves cannot be taken into account and shock 
trajectories are assumed to be radial. Recent three-dimensional 
radiative hydrodynamic simulations result in a dynamic and ther- 
mall y bimodal structure of th e non-magnetic chromospheric lay- 
ers ( We demeyer et al.L [20041) . The simulations show that upward 
travelling and interfering waves and shocks lead to the formation 
of a network of hot upwelling filaments and enclosed cool re- 
gions. The non-magnetic numerical modelling corroborates ob- 
servational findings that the internetwork atmosphere and its be- 
haviour are predominantly determined by the chromospheric re- 
sponse to acoustic oscillations. 



Magnetic fields play 

work grain genera tion .. . - - r 

Remling et al, 1996; iHofmann et al.L 11996 : Steffens e t al 
1996:.Hoekzema & RuttenLll998l:lLites et al.l.ll999i:.De Pontieu . 



a minor role in the in ternet 
von Uexkuell & Kneed. 11995 
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Fig. 1. (a): Theoretical tran smission curve of the Lyot filter (dot-dashed) and the mean spectrum of the Can K line taken from a 
photospheric spectral atlas (iNeckell [19991) . The shaded area indicates the 60 pm FWHM range, (b): Best filtergram taken in the 
Can K line. One minor tick corresponds to 5 arcsec. The filtergram is displayed on a logarithmic scale, (c)-(g): Mean (large) and 
individual (small) magnetograms taken from MDI. One minor tick corresponds to 5 arcsec. 



I2QQ2I: iHoekzema et al.L l2Q02h and are only occasionally involved 
in fo rm of the fla shers ([Bran dt et al.L 119921: iNindos & Zirinl 
119981: iKriiger et al.L I2QQ1I: Ide Wiin et all l2QQ5h . However, mag- 
netic fields appear to play a major role in the higher atmosphere 
where plasma /3 (ratio of gas pressure over magnetic pressure) 
becomes small (iMcIntosh et al.l l200ll) and wave reflection and 
conve rsion between acoustic and magneto-acousti c modes takes 
place ([Rosenthal et al.Ll2002tlBogdan et al.Ll2003h . 

In the present work we investigate the Ca n K chromosphere 
and its distinct fine structure based on an analysis of a long 
time sequence (~6h) of narrow-band (60 pm FWHM) filter- 
grams featuring a spatial resolution of 0.7 arcsec. We concen- 
trate on morphological, dynamical, and evolutionary diff'erences 
between the network and the internetwork. We employ image 
segmentation techniques to establish statistics of properties like 
size and nearest-neighbour distance that can be used to validate 
current and future simulations. 

The paper is organised as follows. The observations are in- 
troduced and the data reduction process is reviewed in Sect. [2l 
and Sect.O respectively. In Sect. IH we go into image segmenta- 
tion processing and describe the applied techniques (a) to distin- 
guish on large- scales between the network and the internetwork, 
and (b) to identify and isolate small-scale chromospheric bright- 
enings. In Sect. [5] the temporal evolution is studied as a result 
from a correlation analysis Section [6l is dedicated to a morpho- 
logical analysis of chromospheric brightenings. The results of 
our investigation are concluded in Sect. [71 

2. Observations 

The observations discussed here have been carried out at the 
Vacuum Tower Telescope of the Kiepenheuer-Institut at the 
Observatorio del Teide on Tenerife, May 7, 1999. A Halle Lyot 
filter was used to obtain narrow-band filtergrams in the Ca n K 
line at 393.3 nm. The width of the filter was set to 60pm and the 
transmission band was centred ~0.017nm away from line core 
on the K2v emission reversal. A ~6 h sequence was recorded with 
a cadence of 6.04 s on a blue- sensitive 1024x1 024 pixel CCD 
camera operated in a 512x512 summing mode. The resulting 
image scale on the detector was 0.35x0.35 arcsec^. The expo- 
sure time was 250 msec. For ima ge stabilisation the VTT cor- 
relation tracker system was used (iSchmidt & Kentischeii[T995l: 
iBallesteros et al.L [l996h . Figure [T] (a) shows the atlas profile of 
the central part of the Can K line taken from Neckel (1999) and 
the theoretical transmission function of the Lyot filter. Figure 
[T] (b) shows the best Ca n K filtergram of the observed quiet 




7 - 



100 



200 
time [min] 



300 



Fig. 2. Rms-intensity fluctuation as a function of time during the 
observation (small dots). Solid: smoothed version of the fluctu- 
ation used to determine all points that are >7.5 % from the local 
average value (large dots). The vertical dashed line indicates the 
time where the best filtergram displayed in Fig. [T] was taken. 



sun region. The observed intensities in the filtergram are dom- 
inated by the emission feature and contaminations from K3 
and Kiv are strongly suppressed due to the rather sharp decline 
of the transmission function. Even in case the K2v reversal is 
absent, the main contribution to the integrated intensity comes 
from the range of heights that contribute to reversal features, al- 
though at a much reduced source function. The influence from 
red-shifted Kiv or blue-shifted K3 is expected to be small within 
the 60 pm filter passband in our observations. To estimate the 
formation height range more quantitatively we performed a pre- 
liminary analysis employing response functions to examine the 
influence of temperature cha nges in the solar atmosphere on the 
filter integrated signal (see iFossum & CarlssonL l2005l) . From 
the response functions computed through a hydrostatic model 
of the solar atmosphere we conclude that our filtergram observa- 
tions basically reflect conditions prevailing in the height range 
350-650 km. Experiments show that this height range is very 
similar in one-dimensional snapshots from simulations of chro- 
mospheric dynamics (e.g. Carlsson & Stein, 1997). For a more 
thorough study of the influence of temperature changes in the 
solar atmo sphere on the resp onse of the filter integrated signal 
we refer to lUitenbroekl (l2006l) . 

MDI on-board SoHO acquired magnetograms and contin- 
uum intensity data in the full-disk mode (2 arcsec spatial res- 
olution) with a cadence of 96 min. These data were averaged 
in time, rescaled and co-aligned to the mean Ca n K filtergram 
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Fig. 3. Left: Mean intensity map (averaged over the whole sequence). Middle: histogram of the mean intensity map (dotted), double- 
Gaussian fit to the data points (full) and single-Gaussian profiles (dashed, dashed-dotted). The shaded area indicates the range of 
relative intensities that is excluded from the statistical analysis. Right: Masks used for the definition of the internetwork (grey) and 
the network (white). Black areas define regions which do neither belong to the network nor to the internetwork according to our 
definition. 



averaged over the whole ~6h observing period (Fig. [3l left). 
Figure [T] (c) shows the resulting mean magnetogram and the 
four frames (d-g) taken during our observing period. The mag- 
netograms show several magnetic elements up to the size of 
-lOarcsec. These structures are persistent during the whole ob- 
servation and coincide spatially with bright regions in Ca ii K 
outlining the chromospheric network (see Fig.[lJ)). The smallest 
structures visible in our data in the internetwork regions span 2- 
3 pixels corresponding to 0.7-1 .0 arcsec. This puts an upper limit 
on the spatial resolution. 

3. Data reduction 

All data were corrected for dark current and off'set and for inten- 
sity inhomogeneities caused by the CCD or any optical surface. 
The guiding was not perfect throughout the whole sequence; 
occasionally, the correlation tracker, working on normal gran- 
ulation, lost its tracking position in moments of bad seeing. 
Therefore the images of the series were aligned by shifting and 
minimising the root-mean- square diff'erence with a precision of 
about 0.2 pixel. To this end, every 100th image was selected as a 
reference image, and these images were aligned first. If the cho- 
sen reference image had poor contrast, we used a neighbouring, 
good-quality image instead. In a second step, the full sequence 
was then aligned with respect to the nearest reference image. 
The post-facto alignment procedure reduced the common field of 
view to about 150x1 50 arcsec^. The remaining image displace- 
ments are below 1 pixel, and are mostly due to distortion within 
the field of view. They do not aff'ect our data analysis. 

To characterise the quality of our data set we computed the 
root-mean- square (rms) intensity contrast for each image, which 
is displayed in Fig. [J] The rms values fluctuate around a mean of 
~9 %, with a maximum and minimum of almost 12 % and 6 %, 
reflecting the best and the worst moments, respectively. 

4. Image segmentation 

This section introduces the image segmentation methods we em- 
ploy to determine statistical properties of internetwork and net- 
work fine structure (see Sect. [6]). To distinguish between inter- 
network and network areas on large spatial scales we employ a 



threshold technique, and to identify and isolate individual bright 
structures on relatively small spatial scales, we apply morpho- 
logical image processing techniques. 

4.1. Definition of network and internetwork 

We use the very diff'erent temporal behaviour of the network and 
the internetwork to distinguish between the two. After normal- 
ising the intensities of each individual filtergram to the mean 
intensity of the image /mean the mean intensity map (Fig. O 
left) was calculated by averaging over all 3483 single filter- 
grams. Motivated by the visual impression that the intensity 
map is dominated by two diff'erent populations of intensity 
we fitted the histogram (Fig. [3l middle) of the mean inten- 
sity map by a double-Gaussian function. The two maxima of 
the Gaussian components, located at 0.97 and 1.02, were used 
as intensity thresholds: pixels darker than < 0.97 and brighter 
than > 1.02 are c onsidered as internetwork and network, re- 
spectively (see e.g. iKneer & von Uexkulli 119931: ISteff'ens et al.L 
119961: iKriiger et al.l 1200 ll) . The right panel of Fig. [3] shows the 
binary mask that results from thresholding the mean intensity 
map with the values provided by the double-Gaussian fit. White 
areas define network regions and grey areas define internetwork 
regions. Black areas are indefinite and are not included in the 
analysis described in Sect. [5] and Sect.[6l If not mentioned other- 
wise in the following all intensities are given in units of /mean- 

4.2. Morphological image processing 

We pursue an image segmentation (or tessellation) approach to 
identify and isolate individual chromospheric brightenings in the 
network and the internetwork. We use an IDL implementation 
of the watershed transformation to segment each filtergram into 
watershed regions (cells) and their boundaries. As a result each 
individual cell harbours one maximum and thus isolates and de- 
fines a single structure. The transformation makes use of the 
morphological watershed operator labelling each watershed re- 
gion with an unique index, and boundaries set to zero. For fur- 
ther details see Dougherty (1992). Prior to the image segmen- 
tation process an optimum filter is applied to each filtergram to 
suppress high-frequency noise and to avoid over- segmentation. 
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Fig. 4. Left: Segmented image as a result from the watershed 
transformation. Right: Binary mask indicating the intensity peak 
position of identified fine structure from the watershed transfor- 
mation. For better visibiHty the positions are marked using more 
than one pixel. 



Figure |4] (left) visualises the tessellation for a sample filter- 
gram. We estimated the background intensity of each identified 
structure inside the cell by averaging over the cell boundary con- 
tours. The number of cell pixels define the cell size and the av- 
erage over these pixels defines the mean cell intensity. From the 
mean background intensity and the peak intensity we derive the 
full-width-half-maximu m size (eff'ective diameter) of the corre- 
sponding structure (see iTritschler & Schmidt. l2002h . Finally, 
each of the identified maxima are connected to its neighbour- 
ing maxima via a mathematical triangulation which returns an 
adjacency list from which nearest-neighbour distances can be 
computed directly. 

In addition we determine all local maxima along each inten- 
sity row and column via numerical diff'erentiation. Local maxima 
in the image are then found by identifying all the pixels where 
local maxima in row and column coincide. A comparison shows 
that these local maxima match very well with those found us- 
ing the watershed transformation. Although the peak intensity is 
determined by the image segmentation process we do not eval- 
uate on the statistics of those but use it only for thresholding 
purposes. 



5. Temporal evolution of the chromosphere 

In the following, we analyse the temporal evolution of the chro- 
mospheric network and internetwork regions. First, we deduce 
characteristic time scales for the network and internetwork, sep- 
arately. Second, the results from the segmentation processing are 
used to visualise and study the appearance of the fine structure. 
Third, we examine the temporal and spatial frequency behaviour 
using the k - v-diagram as well as power maps as a diagnostic 
tool. 



5.1. Autocorrelation Analysis 

In order to specify the temporal evolution in more detail char- 
acteristic time scales are derived from an autocorrelation analy- 
sis. To compare with the numerical results of Wedemever et al] 
(|2004 the whole field-of- view is divided into 22x22 pixel sub- 
fields with 50 % overlap so that the subfield size corresponds to 
the simulation size which covers a -5.6x5.6 Mm^ region of the 
quiet Sun horizontally. Each individual subfield is allocated to 
the network or to the internetwork if 100 % of the pixels coin- 
cide with network or internetwork pixels as defined by a mask 
similar to the one described in Sect. 14.11 but using an intensity 
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Fig. 5. Autocorrelation versus time. Dots: individual autocorre- 
lation functions for all subfields that belong either to the inter- 
network (upper panel, 73 subfields corresponding to the filled 
circles in Fig.[6l right) or to the network (lower panel, 10 sub- 
fields corresponding to the open circles in Fig. [6l right). Thick 
solid: average over all network/internetwork subfields. Thick 
dash-dotted: individual autocorrelation functions. 
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Fig. 6. Left: Evolutionary time scale versus mean intensity. 
Vertical dotted lines indicate the intensity threshold chosen to 
diff'er between internetwork (filled circles) and network (open 
circles) subfields. Diamond: internetwork subfield with compa- 
rably long time scale. Triangles: network subfields with com- 
parably short time scales. Crosses: network subfields with time 
scales longer than the duration of the used data set (see text). 
Right: Evolutionary time scale per subfield. 



map averaged over -75 min (instead of the full -350 min). For 
computational reasons we performed the autocorrelation analy- 
sis for all subfields on the first 750 images (~75 min) only which 
provides a statistical relevance of the autocorrelation function up 
to 25 min. 

In Fig. [5] each single autocorrelation function is plotted as a 
function of time for all internetwork (upper panel) and network 
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(lower panel) subfields, revealing a diversity in amplitude and 
phase. We highlight the autocorrelation of two individual sub- 
fields in the internetwork and network (dash-dotted) to demon- 
strate that their behaviour stays quasi-periodic for a time lag of 
at least 25 min. For both the network and the internetwork, the 
first cycle in the individual and the average autocorrelation peaks 
at ~3min. The following cycles tend to appear with a 3-4 min 
period for internetwork subfields and with a ~5 min period for 
network subfields. 

As an evolutionary time scale we take the time after which 
the autocorrelation decreased to a value of 1 je. For those sub- 
fields where the autocorrelation function does not fall below a 
value of 1 le, the time scale is set to the duration of the chosen 
data set of 4530 sec. Figure [6] shows the derived time scale ver- 
sus the mean intensity per subfield (left), and a spatial map of 
all deduced time scales (right). The time scales show a clear di- 
chotomy between network and internetwork subfields. The in- 
ternetwork time scales group together in the range 43-91 sec 
with a mean of 52 sec. In taking this mean we neglect one out 
of the 74 internetwork subfields that shows a much larger time 
scale of 716 s, most likely due to contributions of unresolved 
magnetic elements. We cannot evaluate the proper mean time 
scales for the network subfields because of the time constraint 
of ~75min (see above). However, network time scales are well 
known to vary subs tantially in the ran ge fro m min to many 
> 1 994: MuUeR [1994: Berger & Titld. 119961: 
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hours (e.g. iLiu et a 
iBerger et al.[ll998l: 



Nisenson et al.L 



Comparison with simulations. Evolutionary time scales based 
on horizontal cross-sections of gas temperature at fixed gt- 
ometrical heights in th e 3D hydrodynamical simulations of 
IWedemeyer et al.l (|2Q04|) of the quiet sun are height dependent 
with values in the range 20-25 sec at chromospheric heights 

i ~ 1000 km) and about 70 sec in the middle photosphere 
WedemeyeR l2003l) . The temperature amplitudes of the model 
chromosphere might be somewhat uncertain due to the simpli- 
fied treatment (grey, LTE), whereas the dynamics and the related 
timescales are connected to already realistically modelled lower 
layers. An indirect qualitative comparison can be done if one 
considers the following eff'ects: (a) the observed intensity origi- 
nates from an extended formation height range instead of a fixed 
geometrical height, (b) intensities refer to corrugated surfaces of 
optical depth instead of plane horizontal cuts like used for the 
model, (c) observations are subject to seeing (spatial smearing) 
and instrumental eff'ects, in particular spatial resolution, (d) the 
simulations cover horizontally a limited range of ~6 Mm only 
and do not account for dynamic phenomena larger than these 
scales, (e) the simulations are purely hydrodynamic. The mean 
value of 52 sec derived from our observations is too short for be- 
ing due to reversed granulation in the middle photosphere (70 sec 
in the model) but, considering the mentioned eff'ects, is actually 
in a plausible range expected for (low) chromospheric layers. 



5.2. Image segmentation analysis 

We turn now to the question whether the internetwork bright- 
enings occur at preferred locations. We use the image segmen- 
tation processing that provides a list of positions of all identi- 
fied brightenings, and compute the cumulative bright point num- 
ber in dependence of integration time. Figure [7] (upper panels) 
shows the result for an integration time of 350 min. We display 
two maps: if all the brightenings are taken into account (left) 
and only those that have intensities > 1.2 /mean (right). For better 
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Fig. 7. Upper panels: Cumulative number of bright points after 
an integration time of 350 min (left) and with a threshold inten- 
sity of > 1 .2 /mean (right). One minor tick corresponds to 5 arcsec. 
Lower panels: Zoom-in (by a factor of two) into the cumulative 
number of bright points (at coordinates X = 48 - 93 arcsec and 
Y = 4 - 49 arcsec as indicated by the white square in the upper 
panel, left) after 3 min (left) and 10 min (right) of integration. 
One minor tick corresponds to 2 arcsec. 



visibility of the internetwork regions the cumulative maps are 
displayed on a logarithmic scale. The top left panel of Fig. [7] 
shows locations that are void of brightenings during the whole 
observing time. These are locations close-by the network and 
around strong magnetic elements where almost no brightenings 
appear nor move to. The map with an intensity threshold value 
of 1 .2 /mean) (Fig.IH uppcr panel, right), shows that there are also 
internetwork regions with a deficit in the occurrence of bright- 
enings like the one in the lower left part of our FOV (extending 
70 arcsec in x-direction and 40 arcsec in y-direction). We cannot 
exclude completely that these internetwork regions directly sur- 
rounding the network are a byproduct of the image segmentation 
process. Network bright points appear to be diffuse in intensity 
because of scattering eff'ects in the solar atmosphere. As a conse- 
quence localised peaks produced by internetwork grains moving 
towards or emerging close by the network are harder to detect 
against the diff'use halo around the network bright points which 
results in a deficit of detected brightenings. The lower panel of 
Fig. [7] shows a zoom-in into the cumulative maps for integration 
times of 3 min (left) and 10 min (right). We recognise a regular 
mesh-like pattern in the internetwork which is still clearly visi- 
ble after an integration time of up to -10 min (a typical granular 
life time). The brigthenings though appear to be mostly clus- 
tered on the nodes of this mesh. Even after an integration time 
of 60 min (not shown) the mesh is still distinctive. On average, 
the cell size of the mesh is ~4-5 arcsec. Since we mark the posi- 
tion of each brightening detected in every successive filtergram, 
the cumulative maps do reflect not only the appearance of a new 
internetwork grain but also the movement of those grains (due 
to seeing and proper motion) that have been detected before. A 
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Fig. 8. Power maps with 3 min (left) and 5 min (right) 
periodicity computed over the full 6 hours. The contour 
lines in the 3 min power map outline the network as 
determined from the intensity map averaged over the 
whole time sequence. The images are scaled indepen- 
dently. One minor tick mark corresponds to 5 arcsec. 
Circles indicate locations with increased signal in the 
5 min band that are co-located with enhanced magnetic 
signal in the average MDI magnetogram (see Fig.[TJ c). 



movie made of 3 min cumulative maps, shows evidence that the 
mesh pattern itself develops and migrates. 

What is the physical mechanism that causes the bright grains 
to appear on an organised mesh? Simulations of solar convec- 
tion show the organisation of strong downdraft fingers on a sim- 
ilar spatial scale (iStein & Nordlundlll989|) . Stro ng turbulence is 
associated with these downdrafts. According to iRimmele et alJ 
(Il995h this turbulence gives rise to acoustic noise, which is a 
potential source for the Can K2v excitation. This would ex- 
plain that the grains appear preferentially on a meso scale dic- 
tated by the downdraft fingers. The locations with a high num- 
ber of accumulated bright points (brightest areas in Fig. [71 
lower left panel) may well correspond to locations of such 
strong downdrafts. The bright points that appear in other loca- 
tions may be generated by weaker turbulence in the more pro- 
nounced intergranular lanes that outline the meso- scale down- 
flow pattern. The latter could explain the slow migration of 
the pattern on time scales corresponding to granular evolu- 
tion while the changes in the location of the more persistent 
bright points corresponds to the times over which the meso- 
scale downflows evolve. Note, however, that there is not nec- 
essarily a one-to-one correspondence between the location of 
the photosp heric excitation a nd th e resulting occurrence of a 
bright point ( Hoekzema et al. 1998: Ho ekzema & RuttenLll998]: 
lLeenaarts~& Wedemeyer-Bohm. .2005.) because the waves may 
not propagate vertically. 

5.3. Power maps 

We compute the power over the ~6 hours at each pixel and apply 
a bandpass filter to extract the contributions from signal with a 
3 min and 5 min periodicity. The bandpass filters are centred at 
5.5 mHz and 3.33 mHz, with a FWHM of 0.4 mHz each. Large- 
scale intensity variations in time are eliminated by a linear fit to 
the data points. To apodise we apply a one-dimensional Manning 
window function that covers 8 % of the data vector length and 
tapers the edges to zero. 

Figure [8] shows the resulting power maps thus visualising the 
well-established eff'ect that internetwork and network respond in 
a diff'erent way to wave motions most likely excited in the photo- 
sphere |Liu & Sheelev, 1971; Cram & Dame, 1983; Dame et al., 
119841: iLites et all Il993[ 1 1994 lJudge et al.l 120011) . The network 
participates predominantly in long-period oscillations (>5 min) 
while in the internetwork short-period oscillations (~3 min, and 
shorter) dominate the picture. 

From the 3 min map we recover that there are regions that ac- 
cording to our definition clearly belong to the internetwork, but 
show a significant suppression of 3 min power. The most con- 



spicuous area is the slanted one between X = 20 - 70 arcsec 
and Y = 0-35 arcsec (also apparent in Figures [3l left, and 
[71 up per panel, left). These voids have bee n observed before 
(e.g. lJudge et all 20011: iKr ijger et al.L l200ll) and named mag- 
netic shadows by lJudge et al. (2001)). The phenomenon is as- 
cribed to the spreading of the magnetic field with height and 
forming the canopy, shadowing its immediate surroundings and 
thus hampering (or modifying) upward propagat ing waves (see 
also [Rosenthal et al.L I2OO2I: iBogdan et al. I l200l . The shadows 
also appear around strong network bright points with a persistent 
magnetic signal as can be seen in the right panel of Fig. [8] and 
in addition as narrow lanes along the network as indicated by 
the map of the cumulative number of bright points (see Fig. [71 
upper panel, left). Inspecting the radial behaviour of azimuthally 
averaged power spectra from the centre of the indicated network 
bright points (BP1-BP4 in Fig.[8l right) outwards we find the ten- 
dency that the integrated power is sharply decreasing with dis- 
tance giving rise to the shadows around the bright points. Power 
increases again to average levels in the internetwork beyond the 
indicated circles. 



6. Morphology 

The results of image segmentation are visualised in terms of 
occurrence distributions separately for the internetwork and the 
network. Our network and internetwork masks (Sect. 14.11) avoid 
regions that show mean intensities between 0.97 < /mean < 1-02 
by definition. For comparison, we also show the occurrence dis- 
tribution representative for all identified brightenings inside the 
FOV. Note that these occurrence distributions include brighten- 
ings with 0.97 < /mean < 1-02. Each occurrence distribution is 
normalised to the total number of identified brightenings. 

Since we do not track each individual identified brighten- 
ing in time, our sample of brightenings is not statistically in- 
dependent. Therefore, and to be more independent from see- 
ing conditions, only those filtergrams contribute to the distribu- 
tions, that show a rms-contrast of 7.5 % higher than the value of 
a smoothed version of the contrast variation at the same time 
step (see Fig. [2]). This leads to a selection of 360 filtergrams 
with a mean cadence of 57 sec (not equidistant). We count only 
those structures that span more than two pixels in both cell size 
and eff'ective diameter, and that have nearest-neighbour distances 
larger than 3 pixels. This leaves us with a sample of 206056 iden- 
tified structures. 
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Fig. 9. Occurrence distributions of the effective diameter (left), the cell size (middle), and the nearest neighbour distance (right). The 
distributions are normalised with respect to the corresponding number of identified brightenings. For each quantity the distributions 
for the network (grey), the internetwork (black), and all identified structures (dashed) are displayed. 



6.1. Size distributions 

The image segmentation techniques (see Sect. ^ are used to de- 
rive three diff'erent size measures: eff'ective diameters, cell sizes, 
and mean distances to next neighbour elements. Figure [9] dis- 
plays the resulting occurrence distributions. They are single- 
peaked and asymmetric with a steep rise to the peak, followed 
by a gradual fall-off' with size. All size measures peak almost 
at equal values for the network and the internetwork: we find 
typical eff'ective diameters of 1.6arcsec for both internetwork 
and network structures. We partly ascribe this to the limitations 
of the image segmentation process. Characteristic cell sizes are 
almost twice as large as eff'ective diameters and peak around 
2.8arcsec and S.Oarcsec while nearest-neighbour distances are 
typically 4.0arcsec and 4.2arcsec for the internetwork and net- 
work, respectively. Arithmetic mean values are larger then the 
characteristic values since they reffect the asymmetry of the cor- 
responding distribution function. The mean value for the nearest- 
neighbour distance in the internetwork amounts to 5 arcsec quan- 
tifying and confirming the impression we got from the the cu- 
mulative bright point number maps (see Fig. [7]) that reveal the 
existence of a ~4-5 arcsec (judged by eye) mesh-like pattern in 
the internetwork. 

6.2. Discussion 

The simulations of ICarlsson & SteinI (1 19971) do reproduce the 
spectral and temporal behaviour of the Ca ii H line profile quite 
well but are one-dimensional and thus cannot predict morpho- 
logical properties. On the other hand, current three-dimensiona l 
radiative hydrodynamic simulations (IWedemeyer et al.L |2004|) 
have to neglect time-dependent Hydrogen ionisation and non- 
grey NLTE radiative transfer, so that no realistic Can K line 
profiles can be synthesised. 

Seeing conditions, diff'erential image motion, and instrumen- 
tal eff'ects inffuence the spatial resolution and thus aff'ect the im- 
age segmentation algorithm. We resolve structures that are larger 
than 0.7 arcsec. Spatially unresolved conglomerates of small 
structures are interpreted as a large single structure, which leads 
to an underestimation of the frequency of small-scale brighten- 
ings. We argue that this is predominantly the case for the in- 
ternetwork. It is most likely that we predominantly detect the 
brightest structures and are unable to resolve their substructure. 
Underestimation of the occurrence of faint structures on the one 



side and not resolving the substructure of the rare brightest struc- 
tures shifts e.g. the peak of the size distributions to larger scales. 
This might explain the similarities found in the size distribution 
functions of the network and the internetwork. To clarify this 
data sets with better spatial resolution are needed. 

What is the physical mechanism that causes the bright grains 
to appear on an organised mesh (see Sect. 15.21) ? Simulations of 
solar convection show the org anisation of strong downdraft fin- 
gers on a similar spatial scale (IStein & Nordlundll 19891) . Strong 
turbulence is a s sociat ed with these downdrafts. According to 
iRimmele et aD (1 19951) this turbulence gives rise to acoustic 
noise, which is a potential source for the Can K2v excitation. 
This would explain that the grains may appear preferentially on 
a meso scale dictated by the downdraft fingers. The locations 
with a high number of accumulated bright points (brightest ar- 
eas in Fig. [71 lower left panel) may well correspond to loca- 
tions of such strong downdrafts. The bright points that appear 
in other locations may be generated by weaker turbulence in 
the more pronounced intergranular lanes that outline the meso- 
scale downffow pattern. The latter could explain the slow mi- 
gration of the pattern on time scales corresponding to granular 
evolution while the changes in the location of the more per- 
sistent bright points corresponds to the times over which the 
meso-scale downffows evolve. Note, however, that there is not 
necessarily a one-to-one correspondence between the location 
of the photospheric excitation and the re sulting occurrence of a 
bright point dHo ekze ma et al.l.[T9 98: Hoe kzema & Rutted. [19981: 
[Leenaarts & Wedemeyer-Bohml 120051) because the waves may 
not propagate vertically. 

7. Conclusions 

We analysed a long time sequence of narrowband Ca ii K filter- 
grams to investigate the morphological, dynamical, and evolu- 
tionary, properties of the solar chromosphere separately for the 
internetwork and the network. Quantifying these properties is 
important for comparison with current and future numerical sim- 
ulations. 

From an autocorrelation analysis on small subfields we find 
average evolutionary time scales of 52 sec in the internetwork. 
Individual values can be found in the range 40-90 sec. The pre- 
dicted time scales of temperature variations in the numerical 
simulations compare favourably with the observed time scales 
of intensity variations in the internetwork notwithstanding the 
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severe simplifications in the simulations. Diff'erences between 
the two time scales may occur because of intensity variations 
are the result of variations of temperatures over a range of op- 
tical depths and these vary with height themselves. In addition, 
seeing and instrumentation aff'ect the spatial resolution leading 
to longer time scales. 

To identify and isolate small-scale features we apply im- 
age segmentation techniques and find evidence for internetwork 
grain patterning on spatial scales of -4-5 arcsec (see Sect. 15.21 
Fig. [7]). We argue that this pattern is related to photospheric 
downdrafts, which occur on similar spatial scales, through the 
enhanced excitation of acoustic waves that these downdrafts pro- 
vide. 

We separated network from the internetwork by a double- 
Gaussian decomposition of the intensity distribution averaged 
over the full duration of the observation. The validity of this 
decomposition is confirmed by the 3 min power map (Fig. [H 
left) in which the power deficit on the network coincides almost 
perfectly with the network mask. There is, however, one region 
(lower left) that appears dark in both the mean intensity map 
(Fig. [H left) and the 3 min power map. We speculate that this 
region contains magnetic field that has not penetrated yet into 
the chromosphere. Either the magnetic field hinders the propa- 
gation of acoustic waves or these waves suff'er mode conversion 
at low altitude before giving rise to Can K emission. The latter 
is also the mechanism proposed for the small-scale shadows we 
see around strong network bright points (Fig.[8l right). 

From a statistical analysis based on results from image seg- 
mentation we find that internetwork grains and network bright 
points appear to have very similar typical sizes (-1.6 arcsec). 
This is somewhat surprising because of their diff'erent physical 
nature as is evident from diff'erences in the power spectra and the 
disparity in their temporal behaviour. We surmise that this has to 
be attributed in part to insufficient spatial resolution. 

For better understanding of the statistical properties in the 
Can K filtergrams and how they are related to physical prop- 
erties in the solar atmosphere both simulations that allow real- 
istic modelling the formation of Can K, and observations with 
higher spatial and spectral resolution are required, preferably in 
the form of two-dimensional spectroscopy. 
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